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Abstract 

We study cosmological dynamics and late-time evolution of an extended induced grav- 
ity braneworld scenario. In this scenario, curvature effects are taken into account via the 
Gauss-Bonnet term in the bulk action and there is also a Chaplygin gas component on 
the brane. We show that this model mimics an effective phantom behavior in a relatively 
wider range of redshifts than previously formulated models. It also provides a natural 
framework for smooth crossing of the phantom-divide line due to presence of the Chap- 
lygin gas component on the brane. We confront the model with observational data from 
type la Supernovae, Cosmic Microwave Background and Baryon Acoustic Oscillations to 
constraint the model parameters space. 
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1 Introduction 



The outcome of supernovae redshift-luminosity distance and also other observational probes, 
that the universe is undergoing an accelerated phase of expansion, has stimulated a lot of 
attempt to explain this unexpected feature. Since dynamics of the universe is described by 
the Friedmann equation which follows from the Einstein field equations in four dimensions, all 
modifications of the Friedmann equation ultimately affect the Einstein field equations too. In 
this respect, modifications to geometric part of the field equations imply some sort of alternative 
geometries, while modifications to matter sector of the theory involve new forms of energy 
densities that have not been observed yet. In fact, the geometric part of the Einstein's field 
equations can be modified to incorporate dark geometry as G^ u + G^ ark ^ = 8nGTj^f\ In 
the same way and within the second viewpoint, Einstein's field equations can be written as 
Gfiu = 8kG(T^ + T^ arfc) ) where and T^ arU) are energy-momentum tensor of ordinary 

matter and dark energy respectively. 

Within the first viewpoint, a well-studied model of modified gravity is the Dvali- Gabadadze- 
Porrati (DGP) braneworld scenario (G. Dvali, G. Gabadadze & M. Porrati 2000; G. Dvali et 
al. 2000; A. Lue 2006 ), in which our four-dimensional world is a FRW 3-brane embedded 
in a five- dimensional Minkowski bulk. The model is characterized by a cross-over length scale 
r c such that gravity is a four-dimensional theory at scales r ^ r c where matter behaves as 
pressureless dust. In the self-accelerating DGP branch gravity leaks out into the bulk at scales 
r r c and the cosmology approaches the behavior of a cosmological constant (C. Deffayet 
2001; C. Deffayet, G. Dvali & G. Gabadadze 2002 ). In the self-decelerating, normal DGP 
branch, gravity leaks in from the bulk at scales r ^> r c , leading to a cosmology which is in 
contrast to the observed late-time acceleration. 

Cosmological dynamics in a braneworld setup that treats Gauss-Bonnet curvature effect 
and the DGP induced gravity in a unified manner, has been studied recently (G. Kofinas, R. 
Maartens & E. Papantonopoulos 2003; R. A. Brown et al. 2005; R. -G. Cai, H. -S. Zhang & A. 
Wang 2005; R. A. Brown 2007; H. Maeda, V. Sahni & Y. Shtanov 2007; J. -H. He, B. Wang & 
E. Papantonopoulos 2007; K. Nozari & B. Fazlpour 2008; K. Nozari & N. Rashidi 2009 ). This 
model which is called GBIG-gravity, is a generalized braneworld scenario that contains both 
UV (ultra-violet) and IR (infra-red) limits in a unified manner: It contains stringy effect via the 
Gauss-Bonnet (GB) term in the bulk action as the UV sector of the theory and Induced Gravity 
(IG) effect which becomes important in the IR limit. The cosmological dynamics and possible 
realization of the phantom-like behavior in this setup are studied recently (M. Bouhmadi-Lopez 
& P. V. Moniz (2008, 2009); K. Nozari and N. Rashidi 2009; K. Nozari, T. Azizi & M. R. Setare 
2009; K. Nozari & N. Aliopur 2009). 

Within the second viewpoint and focusing on the matter sector of the Einstein field equa- 
tions, a well-studied model introduces into a dark energy component called the Chaplygin 
gas (A. Dev, J. S. Alcaniz & D. Jain 2003; L. Amendola et al. 2003; O. Bertolami et al. 2004; 
M. Biesiada, W. Godlowski & M. Szydlowski 2005; X. Zhang, F. Wu & J. Zhang 2006; H. 
Zhang & Z. -H. Zhu 2006; Heydari-Fard & H. R. Sepangi 2008; H. Zhang, Z. -H. Zhu & L. 
Yang 2009; M. H. Mohseni Sadjadi 2009). This model is similar to the DGP model in the 
sense that it is also characterized by a cross-over length scale below which the gas behaves as 
pressureless dust, and above which it approaches the behavior of a cosmological constant. This 
length scale is expected to be of the same order of magnitude as the r c scale of the DGP model 
(M. Roos 2007). Accelerating Chaplygin gas combined with the decelerating braneworld DGP 
model can produce an overall accelerated expansion of the order of magnitude seen (M. Roos 
(2007, 2008a,b); M. Bouhmadi-Lopez & R. Lazkoz 2007). However, both the self-accelerating 
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DGP model in flat space and the standard Chaplygin gas model have problems in fitting with 
present supernovae data (T. M. Davis 2007). 

Recently, astronomical observations with WMAP7 have indicated that the equation of state 
parameter of dark energy can be less than —1 and even can display a transient behavior (E. 
Komatsu et al. 2010 ). A simple way to explain this phenomenon is to consider a non-canonical 
phantom dark energy (R. R. Caldwell 2002), that introduces new theoretical facilities and 
challenges in this field. Phantom fields are a sort of scalar fields with negative sign for the 
kinetic energy term. Indeed, phantom fields suffer from instabilities due to violation of the null 
energy condition, and a phantom universe eventually ends up with a Big Rip singularity (R. 
R. Caldwell, M. Kamionkowski & N. N. Weinberg). Thus it follows immediately that there 
must be some alternative approaches to realize a phantom-like behavior without introducing 
any phantom field in the model. With phantom-like behavior, we mean the growth of the 
effective dark energy density with cosmic time and in the same time, the effective equation 
of state parameter should stay always less than —1. In this regard, it has been shown that 
the normal, non-self-accelerating branch of the DGP scenario has the potential to explain the 
phantom-like behavior without introducing any phantom fields on the brane (V. Sahni & Y. 
Shtanov 2003; A. Lue & G. D. Starkman 2004; R. Maartens & E. Majerotto 2006). The main 
problem associated to this model and its extended versions (L. P. Chimento, R. Lazkoz, R. 
Maartens 2006; M. Bouhmadi-Lopez & P. V. Moniz 2008) is that the effective phantom picture 
always breaks down at some point as the effective dark energy density will cease to be positive 
at a redshift in the past. 

With these preliminaries, in this paper we study cosmological dynamics and late-time evo- 
lution of an extended induced gravity braneworld scenario. In this scenario, curvature effects 
are taken into account via the Gauss-Bonnet term in the bulk action and there is also a Chap- 
lygin gas component on the brane. Our motivation for incorporation of the GB correction 
in the DGP setup is that the DGP model gives just the infra-red modification of the general 
relativity. It would be expected that a consistent DGP braneworld model would have also 
ultraviolet modifications as well, associated to high-energy stringy effects at earlier times. This 
expectation can be fulfilled via incorporation of the GB term in the bulk action. By including a 
Chaplygin gas fluid on the brane, there will be no breakdown of the phantom-like prescription 
in this model. The presence of the Chaplygin gas component also provides a smooth crossing 
of the cosmological constant line in this scenario. 

The paper is organized as follows: In section 2 we start from the action of the scenario 
and derive the corresponding Friedmann equation. Then, we study cosmological dynamics of 
this generalized braneworld setup and investigate possible realization of the phantom mimicry 
without introducing a phantom field. In section 3, we investigate possible crossing of the 
phantom divide line by the effective equation of state parameter and also we study dynamical 
screening of the brane cosmological constant. We constraint the model parameter space by using 
the observational data from type la Supernovae, Cosmic Microwave Background and Baryon 
Acoustic Oscillations data in section 4. Finally, our summery and conclusions are presented in 
section 5. 
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2 The setup 



We start with the following action 



S = ^f d 5 x^^[R^ - 2A 5 + (3[(R^Y - AR^R^ + B^tf****] } 



2k? 




(1) 



v 5 Jy=0 

where the first term shows the usual Einstein-Hilbert action in 5D bulk with the Gauss-Bonnet 
term and /3 is the Gauss-Bonnet coupling constant. The assumption that j3 is no n- negative 
is motivated by string theory, where typically (5 oc ^ 2 string (G. Kofinas, R. Maartens & E. 
Papantonopoulos 2003; R. A. Brown 2007 ). r c is the DGP crossover scale, A is tension of 
the brane, and C m is the matter field Lagrangian on the brane. In our case, C m containes 
ordinary matter (CDM) and a generalized Chaplygin gas component with equation of state of 
the form p c h — — -4r- The cosmo logical dynamics on the brane obeys the following generalized 

Pch 

Friedmann equation (M. Bouhmadi-Lopez & P. V. Moniz 2008, 2009; K. Nozari and N. Rashidi 
2009 ) 



-> 2 4+§)f( 



H 2 



* K 
<£ + - 



( i/2 + |)-|(" + "- + A ) 



(2) 



The bulk contains a black hole mass and a cosmological constant, so that $ is defined as 
$ + 2/3$ 2 = + Jr. The bulk black hole mass originates on the bulk Weyl tensor so that 
the bulk reduces to Schwarzschild-AdS 5 if f3 — and to AdS 5 if (3 = = T. In which follows, 
we restrict ourselves to the case that the bulk black hole mass vanishes, T = , and therefore 



$ + 2/3^ 



As. 
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In this case the bulk cosmological constant is given by A 5 



+ 



12/3 



where 



I is the bulk curvature. Assuming A 5 = corresponding to a Minkowski bulk, for a spatially 
flat FRW brane (K — 0), the Friedmann equation would be as follows 



1 2 



[H 2 - $) 



n 1 2 

r c H 2 - -^{p + Pch + y) ■ 



(3) 



With A 5 = , we find $ = or $ = Our forthcoming arguments will be based on the 

choice $ = 0. Now we look at the conservation equation. In our setup, there is no energy 
exchange between bulk and brane. Therefore, total matter/energy budget on the brane is 
conserved by virtue of the Bianchi's identity so that p tot + 3H(p tot + p tot ) = 0, where p tot and 
Ptot are total pressure and energy density on the brane respectively. Restricting further so that 
conservation holds for two matter components separately, we have for the CDM component 



p m + 3Hp m = 



(4) 



which integrates to the usual p m = p mo (l + z) 3 , whereas for the Chaplygin gas component, the 
continuity equation can be written as 



p ch + 3H(p ch + p ch ) = 0. 



(5) 
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Since p ch = — 4- (M. Roos 2007; M. Bouhmadi-Lopez & R. Lazkoz 2007 ), we find 

Pch 

1 — A - J — 

Pch = (pch)o [A s + a3(1+Q ) j 1+a , (6) 

where A s = ^ 1+a . Using the definition of the redshift parameter z as = j^, and setting 
a — 1 for convenience, we find 

Pc/i = (p c .)o[A + (1 - + ^) 3{Q+1) ] ^ (7) 

This expression supports the interest on Chaplygin cosmologies since it reflects the fact that 
the energy density of such fluids interpolates between dust and a cosmological constant (A. 
Y. Kamenshchik, U. Moschella & V. Pasquier 2001; N. Bilic, G. B. Tupper & R. D. Viollier 
2002; M. C. Bento, O. Bertolami & A. A. Sen 2002). After these preliminaries, we are looking 
for the late-time behavior of the normal branch of this chaplygin GBIG scenario. Defining 

the cosmological parameters as Vt m = , fi A = %fp , fi rc = Ar 2 H -z , = ^/3Hq and 

( \ 2 ^ 

^ c ?i = 3jj%* 4 ; the Friedmann equation on the brane, equation (3), can be expressed in a 
dimensionless form as follows 

e 2 (z) = -2Jq^ c e(z)[i + n p E\z)\ + n m (i + ^) 3 + fi cfe [A s + (i - a s )(i + z) 3(a+1) ]^ + n A , (8) 



where -E(-z) = j^. We assume < A s < 1 and 1 + a > 0. With these conditions, it is possible 
to realize a de Sitter phase at late time (M. Bouhmadi-Lopez & R. Lazkoz 2007 ). As an 
important ingredient of the model, the following constraint equation can be obtained from (8) 
by setting z = 



i + 2^(1 + sip) = n m + n ch + n A . (9) 

Note that based on this relation, the region Q m + Q c h + Q A < 1 in the model parameters space 
is physically unacceptable. Taking the time derivative of equation (9), and using the continuity 
equation for matter on the brane, the Hubble rate can be deduced as follows 



k 3 [E(z)(l + zf 



Vt m 4 



3(l-A 3 )n ch (l+z) a 



[A s +(l-A s )(l+z) :te ]"+i 



Hi 2 E(z) + v^T[l + fyE 2 (z)} + 2^T c ripE 2 (z) 



(10) 



With the previous constraints on A s , and assuming that j3 > ( motivated by string theory as 
(3 oc £l tr i n g) , it is obvious that this relation always has a negative value in the whole physically 

admissible parameters space. This feature is plotted in fig l 4 . Since H < 0, the Hubble 
parameter decreases as the brane expands, consequently there is no super-acceleration in this 
braneworld universe. Therefore, the brane does not hit a big rip singularity as its fate. Note 
also that H vanishes when z — > — 1, while H is positive at this limit. This is a reflection of late- 
time de Sitter character of the model. It has been shown that it is possible to have Big Freeze 
singularity in FRW universe filled with a generalized Chaplygin gas (A. V. Yurov et al. 2008; M. 



4 The numerical values of cosmological parameters used to plot figures of this paper are taken from Table 2. 
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Bouhmadi-Lopez, P. F. Gonzalez-Diaz & P. Martin-Moruno 2008; M. Bouhmadi-Lopez et al. 
(2009,2010) ). So, essentially appearance of this type of singularity in our framework is probable 
too. Nevertheless, existence of induced gravity on the brane and the Gauss-Bonnet term in the 
bulk action may help to overcome this difficulty This issue needs further investigations and 
we are going to study its separately. The deceleration parameter, q depends on H through the 




Figure 1: Variation of jp with redshift. The parameter values to plot this figure are taken from 
Table 2. 



relation q = — [-j^ + 1], which can be expressed in the following form 



q 



i- 3 - 



r(l+£l' 
I E(z) 



'-] (pm + 



3(l-A 3 )n ch {l+z) 3a 
[A a +(l-A s )(l+z)^]^Ti 



2 [E{z) + v^T(l + VtpE 2 (z)) + 2^rn p E\z)} 



:n) 



In figure 2, the dimensionless deceleration parameter q is plotted versus the redshift for a 
fixed set of the parameters as are presented in table 2. In this model, the universe enters the 
accelerating phase at z ~ 0.13. 



3 Crossing the phantom divide 

In this section, we show that a phantom-like behavior can be realized on the brane, in a relatively 
wider range of redshifts than previously formulated models ( see for instance R. Maartens & 
E. Majerotto 2006; L. P. Chimento, R. Lazkoz, R. Maartens 2006; M. Bouhmadi-Lopez & P. 
V. Moniz 2008). This phantom-like behavior occurs without including any phantom matter. 
The phantom-like prescription is based on the definition of an effective energy density which is 
corresponding to a balance between the cosmological constant and geometrical effects encoded 
in the Hubble rate evolution (V. Sahni & Y. Shtanov 2003; A. Lue & G. D. Starkman 2004; V. 
Sahni 2005). This behavior is based on the definition of an effective energy density p e ff, which 
increases with cosmic time and an effective equation of state parameter always less than — 1. 
More precisely, the effective description is inspired in writing down the modified Friedmann 
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Figure 2: Variation of q versus the redshift. Parameters values are taken from Table 2. 



equation of the brane as the usual 4D Friedmann equation(A. Lue & G. D. Starkman 2004; M. 
Bouhmadi-Lopez & P. V. Moniz 2008 ), so that 

H 2 = ^(Pm+Peff). (12) 

Using equations (8) and (12) , we find 



Peff 



3H. 



2 v /o>(z)[l + n p E 2 {z)\ + Q ch [A s + (1 - A s )(l + zf^]T^ + fi A l (13) 



By definition, phantom-like prescription breaks down if p e ff < 0. Figure 3 shows variation of 

P^efP versus a and redshift in this model. The phantom-like behavior can be realized for all 

values of a that 1 + a > 0. Figure 4 shows variation of p^J/P versus A s and redshift in this 
model. The phantom-like behavior can be realized if < A s < 1. Based on this analysis, p e // 
is always positive and grows by cosmic expansion, and this is a typical phantom-like behavior. 
Note that in the absence of the Chaplygin gas component on the brane, the effective phantom 
picture breaks down at a redshift in the past where the effective dark energy density becomes 
negative. The effect of the Chaplygin component as a new ingredient added to the GBIG 
braneworld scenario, is that the effective phantom-like picture has no break down in this case. 
The main point here is the fact that it is possible essentially to save phantom-like prescription 
for a wide range of redshifts in this case. Now, the effective equation of state parameter cu e //, 
can be defined using the conservation equation of the effective energy density on the brane 

Peff + 3H(l + u eff )p eff = 0. (14) 

A straightforward calculation shows (M. Roos 2007; M. Bouhmadi-Lopez & R. Lazkoz 2007 ) 

. 9M \(1-4)(1 + Z )^) 

Pch [A 8 + (1 - A 8 )(l + z)W+ a )]^ ' 



Then the effective energy density evolves as follows 



Figure 3: Variation of the effective dark energy versus a and redshift. Parameters values are taken 
from Table 2. 




Figure 4: Variation of the effective dark energy versus A s and redshift. Parameters values are taken 
from Table 2. 
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Peff 



m 2 (3H(i + zf(i-A s )n ch (i + 



K 2 



\3a 



[A s + (l-A s )(l + ^)3«] — 

2y^T J B(z)[l + nfsE 2 (z)}HH + 2^\HH VL l3 E 2 (z) 
Using equations (13) and (16), we find 



(16) 



1 + U ef f 



3H(l+zf(l-A s )n ch {l+z) :i °' 
[A a + (l-A s )(l+z)3^]^TT~ 



+ 2y/n^ c E(z)[l + n p E 2 (z)]HH + 2y/n^ c HH o tt E 2 



3E(z) - 2^TE^)[l + n a E*(z)} + tt ch [A s + (1 - A s )(l + z)3(«+D]— + fi A 



(17) 

Figure 5 shows the plot of l+u e ff versus the redshift. The universe enters to the phantom phase 
smoothly at z ~ 0.12. It is important to note that this setup realizes a smooth transition to the 
phantom phase, the so called phantom- divide line crossing. We note that this smooth crossing 
behavior cannot be realized in the absence of the Chaplygin fluid component. Indeed, without 
this term, the effective equation of state parameter blows up in the past and as mentioned 
previously, the effective phantom description breaks down. 



0.1 
0.05 


-0.05 
-0.1 




0.2 0.4 0.6 0.8 1 



Figure 5: 1 + w e ff versus the redshift. This model realizes (as a result of the existence of the 
Chpalygin component on the brane), a smooth crossing of the phantom divide line. Parameters values 
are taken from Table 2. 



Now we study the phantom-like behavior in another perspective: dynamical screening of 
the brane cosmological constant (A. Lue & G. D. Starkman 2004 ). The normal branch of this 
model can be described by the following Friedmann equation 

H 2 = —(Prn + Pch) + ~ ~ - 1 + -PH 2 ) (18) 

3 3 r c 3 

where A is the brane cosmological constant. We rewrite this equation as follows 

H 2 = ^(p m + Pch) + ^A (e//) , (19) 
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where by definition 



A (e//) = A _ H(l + lf3H*) 
3 r r 



(20) 



This equation means that the brane is extrinsically curved so that shortcuts through the bulk 
allow gravity to screen the effects of the brane cosmological constant at Hubble parameters 
H ~ r^ 1 where r c is the DGP crossover distance. As a distinctive feature of this model, 
curvature effect via the Gauss-Bonnet term and presence of the Chaplygin component on the 
brane contribute in the dynamical screening of the brane cosmological constant. While the role 
played by the GB term is obvious, the role of the Chaplygin matter is hidden in the definition 
of the Hubble parameter. This feature leads to a considerable difference relative to the pure 
DGP case studied in (A. Lue & G. D. Starkman 2004 ). Figure 6 shows the difference between 
two scenarios in this respect for some values of Qp. Alternatively, we can define also 




£2=0 £2 =0.01 n =0.02 

p p p 



Figure 6: A.( e ff) versus the Hubble parameter for different values of 17^ 



„(.//)_ A HQ- + l(3H^ 

Pde — T7 

3 r. 



(21) 



Figure 7 shows variation of versus z in our proposed setup. It is always positive and 

grows with cosmic expansion; a typical phantom-like behavior. 



4 Observational Constraints 



In this section we study the constraints imposed on the GBIG chaplygin model parameters 
using the observational data such as the gold sample of SNIa combined with the information 
from the BAO measurement by SDSS and the CMB shift parameter from WMAP7 observations. 
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0.1 0.2 0.3 0.4 0.5 0.6 



Figure 7: Pde versus the redshift. Parameters values are taken from Table 2. 



A: SNIa data 

In this section we use the 156 new gold sample supernovae la data compiled in (A. G. Riess 
et al. 2007) to fit the model. This observation directly measures the apparent magnitude m 
of a supernovae and its redshift z. The apparent magnitude m is related to the luminosity 
distance of the supernovae through 



H = m — M = log 10 d L + 5 log 



10 



cH 
Mpc 



+ 25, 



(22) 



where M is the absolute magnitude that is believed to be constant for all Type la supernovae 
and dj^ is the luminosity distance that in a flat universe can be expressed as follows 



di 



:i+z) 



* dz' 



E(z'\ i 



where 9 denotes the model parameters. From equation (9), the E(z) is given by 



E(z) 



i 1 - \2Vt R Vt 



Pa 



- 1 



(23) 



(24) 



where 



P = im p tt Tc + 54QOjfi rc - 1 + 6^3fi^f2 rc ^16fi^fi2 c + im p VL r . c Q + 24Q r QlQ 2 - Q rc - Q, 
and 

Q = fi m (l + zf + Q ch [A s + (1 - A s )(l + z) 3 ^]— + Ov. 

For the statistical analysis of the supernovae data we perform the x 2 statistics for the model 
parameter that is calculated as 



1=1 



(25) 
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Table 1: Priors on the parameter space used in the likelihood analysis. 



Parameter 


Prior 




^ch 


0.00 - 1.00 


Top Hat 




0.00 - 1.00 


Top Hat 




0.00 - 1.00 


Top Hat 




0.02265 ±0.00059 


Top Hat (BBN) 



In this relation, N = 156 is the number of SNIa data points, n obs is the observed distance 
modulus and the <7j is the uncertainty in the observed distance modulus, which is assumed to 

be Gaussian and uncorrelated so that the likelihood is proportional to e 2 . The parameter 

M = 51og 10 (^) + 25, 
Mpc 

is a nuisance parameter and should be marginalized (integrated out). This can be done by 
following the techniques described in (M. S. Movahed, M. Farhang & S. Rahvar 2009 ) and can 
be performed by expanding the Xsn °f equation (25) with respect to M as 

xl N {6) =A- 2MB + M 2 C (26) 

where 



N 



i=l °i 



and 

^ ) = ^^)-M^M = 0,g) ; c = j:\. (28) 

i=l a i i=l °i 

Equation (25) has a minimum for M = ^ at Xsn(@) = ^(^) — ~ ''^P ■ Using this equation the 
best fit values of model parameters as the values that minimize XsNia(@) can be obtained. 

For the likelihood analysis we marginalize the likelihood function L = exp(— x 2 /2) over h 
and Q m = Q b where is the energy density of the barionic matter. We adopted Gaussian 
priors such that h = 0.705 from the WMAP7 (K. Komatsu et al. 2009 ) and tt b = 0.0456. 
Table 1 summarizes these priors. 

B: CMB shift parameter 

We use the CMB data from WMAP7 observation that includes the shift parameter TZ and 
the redshift of the decoupling epoch z*. The shift parameter 1Z is related to the angular 
diameter distance to the last scattering surface, the comoving size of the sound horizon at 
and the angular scale of the first acoustic peak in the CMB power spectrum of the temperature 
fluctuations. The CMB shift parameter is approximated by (see for more details R. Lazkoz & 
E. Majerotto (2007 ). ) 

n= ^H r(z*), (29) 
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where r(z) is the co moving distance to the redshift z defined by 



The constraints on a typical model using CMB shift is obtained from minimization of the 
quantity 

2 _ [n obs - n th } 2 

XCMB — ~2 ' 

°CMB 

where !Z b s is the observed value of the CMB shift parameter performed from WMAP7 obser- 
vation (K. Komatsu et al 2009) and 1Z t h is corresponding to the theoretical value calculated 
from equation (29). 

C: BAO observation 

The baryonic acoustic oscillation (BAO) peak detected in the SDSS luminous red Galaxies 
(LRG) is another tool to test the model against observational data. BAO are described in 
terms of a dimensionless parameter 



A 



Hodl(z;e) !V3_ ^ 



H(z sdss ; 6)z 2 sdss {l + z sdss f 



The best fit values of the model parameters can be determined by constructing minimization 
of the quantity 

XsDSS — ~2 • \ 66 ) 



^sdss 



/ \ -0.35 

The observed value A bs from the LRG is A b s = 0.469( ± 0.017 measured at z sdss = 

0.35 (R. Lazkoz & E. Majerotto, 2007 ). Here n s = 0.963 is the spectral index as measured 
by WMAP seven year observations (K. Komatsu et al 2009 ). It is important to note that the 
above observational data are uncorrelated form each other, since they are given by different 
experiments and methods. Then we can construct a joint analysis of them as 

Xtot — XSN + XCMB + XSDSS- l^J 

With these preliminaries, we have obtained the best fit parameters of the GBIG chaplygin model 
for SNIa data, the joint analysis of the SNIa and CMB, and finally the combined analysis 
of the total datasets. Table 2 indicates the results of the observational constraints on the 
free parameters of the model. In figure 8, using the best fit values of the model parameters 
and observational data from gold sample, we have compared the theoretical predictions of 
the distance modulus in our model. Figure 9 shows the marginalized relative likelihood with 
respect to parameter A s fitted with SNIa gold sample, SNIa+CMB and SNIa+CMB+SDSS 
experiments. We plot the joint confidence interval of Q rc and As in figure 10 in lcr, 2a and 3<r 
level of confidence for the mentioned datasets. 

Finally and for completeness of discussions, we compare our Chaplygin GBIG scenario with 
ADGP and ACDM ( see for instance (W. J. Percival et al. 2007 ) for a general framework). 
Using the data set in the previous subsections, we obtain the best fit parameters of the ADGP 
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Table 2: The best values for the parameters of GBIG chaplygin model by fitting with SNIa 
Gold sample, SNIa+CMB and SNIa+CMB+SDSS experiments in a flat background. 



Observation 


f2 rc 




A s 




&>ch 


a 


Xminl Nd.o.f 


SNIa( Gold Sample) 


0.43 


0.11 


0.50t£i£ 


1.45 


0.99 


0.99 


0.923 


SNIa(Gold Sample)+CMB 


0.53 


0.31 


i9 +u - i5 

u - iz -0.09 


1.70 


0.99 


-0.98 


0.943 


SNIa(Gold Sample) +CMB+SDSS 


0.51 


0.41 


0.099t£& 


1.79 


0.99 


-0.98 


0.992 



and ACDM models in table 3. Comparing tables 2 and 3, we see that in our model the best 
fit values ( which are given by x 2 l^d.o.f) are more reasonable than ADGP and ACDM ones on 
observational ground. In figure 11 we have plotted the marginalized relative likelihood with 
respect to parameter Q m fitted with SNIa gold sample, SNIa+CMB and SNIa+CMB+SDSS 
experiments. In summary, as tables 2 and 3 show, our Chaplygin GBIG scenario has better 
agreement with observations than ADGP and ACDM. In other words, existence of a Chaplygin 
component on the DGP brane brings the scenario to be more viable in observational viewpoint 
than the pure ADGP or ACDM. This is a result of wider parameter space available here which 
leads to further degrees of freedom. 



Table 3: The best values for the parameters of ACDM and ADGP models by fitting with SNIa 
Gold sample, SNIa+CMB and SNIa+CMB+SDSS experiments in a flat background. 



Observation 


ACDM model 


ADGP model 


SNIa 


n m = 0.35l™4, X 2 mt JNd.o.f = 0.921 


fl m = 0.007+ - 08 , fl rc = 0.22, xLJNd.o.f = 0.914 


SNIa+CMB 


On = 0.341"$, xLJNd.o.f = 0.926 


n rn = o.33i^, n rc = o.oo, xinJNd.o.f = 0.926 


SNIa+CMB+SDSS 


Q m = 0.32+^, xinJNd.o.f = 0.960 


n m = 0.34+^, n rc = 0.00, ximJNd.o.f = 0.966 
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Figure 8: Distance modulus of the SNIa Gold sample versus the redshift. Solid line shows the best 
fit values with the corresponding parameters of h = 0.70, fi rc = 0.43 and As = 0.50 in la level of 
confidence with Xmin = 141.368 for GBIG chaplygin model. 




Figure 9: Marginalized relative likelihood with respect to parameter As fitted with SNIa gold sample, 
SNIa+CMB and SNIa+CMB+SDSS experiments. The intersection of the curve with the horizontal 
solid line is corresponding to the bound with la level of confidence. 
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Figure 10: Confidence interval of Q rc and As in la, 2a and 3a level of confidence for the SNIa gold 
data (A), SNIa+CMB joint data (B) and SNIa+CMB+SDSS joint data (C). 
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Figure 11: Marginalized relative likelihood with respect to parameter m fitted with SNIa gold sam- 
ple, SNIa+CMB and SNIa+CMB+SDSS experiments for ACDM (left) and ADGP (right) respectively. 
The intersection of the curve with the horizontal solid line is corresponding to the bound with la level 
of confidence. 



5 Summary and Conclusion 

In this paper, we have considered a braneworld induced gravity scenario which contains cur- 
vature effect of the Gauss-Bonnet term in the bulk and a Chaplygin gas component in the 
brane action. We have investigated the cosmological dynamics and the late-time behavior on 
the brane. We have shown the possibility of getting accelerated expansion in this scenario. We 
have studied the effective phantom-like behavior on the brane. This behavior can be realized 
in this scenario and the universe enters into the phantom phase at z ~ 0.12. We note that 
this effective behavior could be deduced also in some braneworld models (such as ADGP and 
pure GBIG scenario) with number of parameters less than our proposed model. However, in 
such models the effective phantom picture breaks down, and the effective equation of state 
parameter blows up at a redshift in the past. In our model, inclusion of the Chaplygin gas 
fluid on the brane saves the phantom-like prescription for a wide range of redshifts. In fact, 
there is no break down of the effective phantom picture in the presence of the Chaplygin com- 
ponent on the brane. From another perspective, existence of the chaplygin gas component on 
the brane, leads to a smooth crossing of the phantom divide line by the effective equation of 
state parameter of the model. We have studied also the notion of dynamical screening of the 
brane cosmological constant in this generalized setup. Finally, We have confronted the model 
with observational data from type la Supernovae, Cosmic Microwave Background and Baryon 
Acoustic Oscillations to constraint the model parameter space. The results of comparison of 
the model parameters space with the observational data are summarized in tables 1 and 2, and 
corresponding figures 8-10. We have shown that up to the analysis on the parameters space 
performed here ( table 2 and 3), this model has better agreement with observations than the 
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ADGP and ACDM. In other words, existence of a Chaplygin component on the DGP brane 
with curvature effect, brings the scenario to be more viable in observational viewpoint than the 
pure ADGP case. This result is due to the wider parameter space and also further degrees of 
freedom accessible in the model discussed here. 
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